The Mass of a Spin Vortex in a Bose-Einstein Condensate 
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In contrast to charge vortices in a superfluid, spin vortices in a ferromagnetic condensate move 
inertially (if the condensate has zero magnetization along an axis). The mass of spin vortices 
depends on the spin-dependent interactions, and can be measured as a part of experiments on how 
spin vortices orbit one another. For Rb*^ in a 1 pm thick trap, m„ ~ 10~^^kg. 



Vortices, with their long-life and concentration of en- 
ergy, often provoke comparison to particles. But does 
their motion fit the analogy? Ordinary vortices in a 
fluid or superfluid do not move inertially, as particles do, 
because their motion is Magnus-force dominated. For 
example, in the absence of a force, they do not move. 
When a force is applied, they move perpendicular to it, 
a situation which is described by first order differential 
equations In fact, the motion of a pair of vortices is a 
miniature version of how Descartes explained the mo- 
tion of planets, with the sun causing the ether to whirl 
around, dragging the planets at the same speed. On the 
other hand, spinor superfluidsj^, 0? ] made out of laser- 
cooled atoms can have spin-current vortices. These vor- 
tices will be argued to obey Newton's laws at low speeds; 
in particular they have a mass, which determines their 
resistance to being accelerated. 

A mass can also be defined for a vortex in a super- 
conductor or a superfluid and it may play a role in de- 
termining the oscillation frequency of vortex lattices and 
the tunneling rate of vortices 0. Observing the mass for 
such vortices is much more subtle than for spin vortices, 
though, because the vortex inertia is nearly overcome by 
the Magnus force. The picture of a moving spin vortex 
described here is inspired partly by the motion of vortex 
rings[7| and Ref. ^]'s study of that motion based on the 
Gross-Pitaevskii equations. 

We will focus on the case of spin 1 ferromagnetic atoms 
in a two-dimensional condensate. A magnetic field is 
applied along z to stabilize the spin vortices, via the 
quadratic Zeeman effect. (We will take the condensate 
to be in the xy plane for definiteness, but the spin and 
spatial coordinates can be chosen independently.) The 
ground states of these atoms have the form 





(1) 



The values ni,no,n_i correspond to the optimal pro- 
portions of the three spin states (the total density is 
n — ni + hq + n^i), and are fixed. If the magnetiza- 
tion along the z-axis, Mz = ni — n_i, vanishes, then the 
spin is in the ^-plane, as is preferred by the quadratic 
Zeeman effect [5| . The overall phase of the wave function 
(which is not observable) is 9c, and —9s is the azimuthal 



angle of the spin. (To see this, calculate < > and 
< Sy > for this state.) The latter angle can be measured 
by scattering polarized light off the condensate. 

There are two types of vortices in a spinor condensate 
in a magnetic field: a charge vortex, described by 9c — 
±(/), = and a spin vortex described by 9c = 0,9s = 
±0 (and observed in a rubidium-87 condensate^). Here, 
(j) is the azimuthal angle centered on the vortex core. Fig. 
[1] illustrates these vortices: 




V's = 




(2) 



The spin texture is uniform around a charge vortex (ex- 
cept near the core); the spin direction rotates by 360° 
clockwise or counterclockwise in a spin vortex. See Fig. 

m 

A spinor condensate is somewhat similar to a mixture 
of several species of atoms. From this perspective, a spin 
vortex is a bound state of two opposite vortices in two 
different components and a charge vortex is a bound state 
of three vortices. 

The velocity fields in the components of a spinor con- 
densate are related: 



ui^uc + us; 1*1= tie; u^i = uc~us; 



(3) 



The flow in the middle component (for spin 1) is the mean 
of the flows in the other two components. As a simple ex- 
ample, the atoms near a charge vortex flow with velocity 
Um — -p-V^m 0- They all move in the same direction 
(see Fig. [T^), meaning that there is a net transport of 
mass. On the other hand, the atoms near a spin vortex 
do not carry any mass if the condensate has zero mag- 
netization Mz = ni — n_i, because the atoms of spin 
±1 move at equal and opposite speeds. However, there 
is a net transport of spin, since the amount of angular 
momentum carried counterclockwise depends on the cur- 
rent of TO = 1 atoms moving counterclockwise and of 
TO = — 1 atoms moving clockwise. Vortices can be clas- 
sified by the amounts 9c and 9s wind by, Qc and Qs 
respectively. (In general, the charge and spin currents 
are given by Js = nuc + MzUs and Js = QzUs + MzUc 
where qz — ip'' Sztp. If Mz = 0, charge vortices have only 
charge currents and spin vortices have just spin currents.) 
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FIG. 1: Charge and Spin Vortices of strength 2n. a), c) The 
currents of the spin states around charge and spin vortices, 
b), d) The spin textures around charge and spin vortices. 

Now the force on a vortex moving at speed v and 
tossed about by a flow of charge and spin is given by 
F = —Hz X [X] "mQm(itm " v)]. Thc tcrm for a given 
value of m describes the lift force on the vortex in that 
component: according to Bernoulli's principle higher 
velocities correspond to lower pressures, so the vortex 
moves to thc side where its velocity field is pointing in 
the same direction as the relative velocity of the fluid. 
For a spin vortex of charge Qs the total force is 



p = Fs ~ Qsv X hM^z, 



(4) 



where Fs is a force produced by the background spin 
current and the second term is the Magnus force, or lift, 
responding to the vortex's own motion. For a charge 
vortex of charge Qc, the force (assuming — 0) is 



Qc'V X hnz, 



(5) 



where Fc is produced by charge current. 

Idiosyncracies of Vortex Motion: Let us compare the 
motion of charge and spin vortices. A charge vortex in 
stationary fluid cannot drift along a straight line, because 
the lift force would push it sideways. Furthermore, the 
equations of motion take the form of first order differen- 
tial equations for charge vortices when the inertial term 
p in Eq. ([4]) is assumed small: 



dr 

dt 



1 
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-z X Fc] 



(6) 



motion is perpendicular to applied forces. The vortex 
velocity required by Eq. ^ turns out to equal the back- 
ground flow speed, in accordance with Descartes's con- 
ception of planetary motion. The motion of charge vor- 
tices is determined once their initial positions are given. 

In contrast, a spin vortex behaves in a Newtonian way, 
as long as the condensate has zero magnetization. In 
the absence of spin current, Eq. ^ implies "Newton's 



first law of spin vortices" : a spin vortex in a charge cur- 
rent can move at any constant speed. There is no lift to 
push the vortex off course in an unmagnetized conden- 
sate because the component vortices ■0 rotate in opposite 
directions (see Fig. [TJ:). Now if there is a spin current 
the lift-forces on the component vortices from the coun- 
terpropagating flows add to produce a nonzero Fs. The 
solution to p = 0, V — , , does not make sense if 
Mz — 0. Therefore the inertial term cannot be neglected 
and "Newton's second law of spin vortices" results: a 
spin-vortex in a spin current must accelerate, at a rate 
proportional to Fs. The spin force has an electrostatic 
form: 
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(7) 



One can introduce the vortex mass my by assuming that 



p = rUyT 



(8) 



at least at low speeds. Hence, the equation of motion 
prescribes the acceleration: 



' dt^ 



(9) 



These Newton's laws describe spin vortices in conden- 
sates of any spin, as long as = o[l9l|. 

Phase Space: The phase space of N vortices seems 
likely to be 2N dimensional, given that the basic equation 
for the evolution of a spin-1 condensate is the following 
first-order differential equation: 



ih 



dip 
dt 



2M 



dV 

dip"! 



(10) 



Here V is the potential energy of the atoms. Only the 
spatial coordinates of the vortices seem essential for pa- 
rameterizing the superfluid wave function. 

For charge vortices, the phase space is 2N dimen- 
sional. The motion is even described by Hamilton's 
equations [loj. with the y coordinate of each vortex acting 
as the momentum conjugate to x\ (There are surprising 
consequence for the thermodynamics of vortices 11, 1^.) 
The precise conjugacy relation reads 



{x,y} = ± 



1 



(11) 



where the sign depends on the direction of circulation. 

However, Newton's second law implies that A'' spin vor- 
tices have a 47V dimensional phase space. For ordinary 
objects, the momenta exist only in an abstract space: a 
photograph of a ball in the air does not reveal its des- 
tination. Thc future of a spin vortex will also depend 
on which way it is moving at a given time, so Eq. pO|) 
implies that something about the wave function which is 
observable in a single photograph can be used to deduce 
the vortex's velocity. 
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Vortex in tp^ 



Vortex in 



FIG. 2: Seeing the momentum of a spin vortex. The phases of 
tp±i increase in the directions of the arrows. The component 
vortices are pulled apart by the lift force to a distance Dy oc 
Px- The stretching is limited by the phase-locking produced 
by the spin-dependent interaction. 



To guess the tell-tale trait, note that the lift forces on 
the component vortices of a moving spin vortex pull them 
in opposite directions, but they are restrained from drift- 
ing apart completely (see below). Thus, the stretching 
of the spin vortex increases with its speed (see Fig. [2]). 
The momentum p can be discerned from a snapshot if it 
is proportional to the stretching: 



p = —hniz X D 



(12) 



where D points between the component vortices. This 
relation is consistent with the canonical mechanics of 
the component vortices[i3]. According to Eq. (fTT|) 
{x±,y±} — where {x±,y±) is the location of the 

component vortex with circulation ±2tt. Now the center 
of mass coordinates x and y of the spin vortices have the 
ordinary Poisson bracket of zero, but x is conjugate to 
Px = nih{y+ - y_): 



{x,Px} = 1- 



(13) 



Eq. is correct for a dipole of vortices in a single- 
component superfluid also, but the velocity of the dipole 
increases as its momentum (the distance between the vor- 
tices) decreases! This unusual fact can be traced back 
to the logarithmic interaction energy E cx InD. Spin 
vortices behave more normally because the confinement 
energy is quadratic, so v = ^ increases linearly with p. 

Eq. (I12p holds when the phase-only approximation 
applies 14]. The general formula for the momentum is 
p — —z X JJ cPu r curl J, where J is the mass current 
(see [lEj). In the phase-only approximation, curl J is a 
dipole of delta- functions, leading to Eq. (fT2|l . 

Confinement and the Vortex Mass: Spinor condensates 
can have special symmetries that ordinary mixtures (e.g., 
of two different atoms) cannot have[l6!|. But the higher 



symmetry is not the only special thing about spinor con- 
densates. Unlike in a mixture, the atomic states can turn 
into one another in a spinor condensate, as long as the 
angular momentum does not change: 1 H — 1 ^0-1-0. 
This "coherent spin-flipping" process is contained in the 
potential energy operator, V = ^a.mi„i2 ■ l^miPlV'mjP : 

+2/33?(^lV'liV'o') + #^52^^ - fiH\^. 

Coherent spin-flip ping locks the phases of the spinor 
components together 17,[l8| and this keeps the two parts 
of a spin vortex bound together. A nonzero phase A = 
{6 1 + — 26*0) costs energy 



Sfiip = 2/3noy/nin_i cos A, 



(14) 



where /3 < is the spin-dependent interaction parame- 
ter. Since the energy cannot depend on a phase unless 
particle numbers are not conserved, it makes sense that 
this comes from the spin- flipping reactions. 

This "coherent chemistry" interpretation of Eq. 
is complemented by an argument which views the atoms 
as classical magnets: the magnetization in the xy-plane 
depends on the relative phases between the spin compo- 
nents and decreases as A increases. This goes against the 
ferromagnetic propensities of the atoms. 

Now the potential energy of the stretched vortex 



E ' 



(15) 



can be reinterpreted as kinetic energy on account of Eq. 
(fT2l) . At low speeds, E = so 



2m„ 



(16) 



The second expression is obtained by relating (3 to the 
width of the condensate w and the scattering-length dif- 
ference Aa flS ]. In this case, ruy is of the same same 
order as the mass of the atoms in the vortex core, nl'^M, 
where Im is the magnetic healing length. For rubidium 
in a 1/im wide trap with Aa ^ lA, ~ lO^'^^ kg. 

Measuring the Mass: Now consider the consequences 
of spin vortices' inertial motion. Spin vortices of opposite 
signs orbit around one another whereas opposite charge 
vortices push each other along parallel lines As for 
planets, the orbits have different shapes depending on 
the initial momenta of the vortices. Specifically, when 
the two vortices move on a circle, attracting each other 
according to a ^ force law, the period is proportional to 
its radius. Balancing the attraction Eq. ([7]) against the 
centrifugal force p = pu, gives 



vp{v) 



M 



(17) 



where v is the speed of both vortices. This remains true 
when p{v) is nonlinear. The left hand side is an increasing 
function of v, so v = Vdrc is determined. This speed is 
independent of the orbit radius. 
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Measuring Vdrc will give an estimate for the vortex 
mass appearing in Eq. ([8]). However, the linear relation 
for p{v) is not reliable at the speed Vdrc because it is on 
the order of the speed of spin-waves in the condensate. 

A more accurate way to measure the mass of spin vor- 
tices is to observe their motion when the magnetization 
is not zero, but is small. Then the lift force due to the 
vortex's motion looks like the Lorentz force from a small 
magnetic field, Ba = —Alzhz. A single vortex will there- 
fore follow cyclotron orbits with the period 

If the magnetization is 5% and the other parameters of 
the condensate are those given above, this period comes 
out to be .3 sec. 

Limiting Velocity: A final idea for an experiment is 
to study the motion of a rapidly moving vortex. Such a 
measurement allows hypothetical inhabitants of the su- 
pcrfluid to measure the velocity of their "ether." Spin 
vortices can move inertially only up to a certain velocity 
relative to the condensate. 

Imagine pushing a spin vortex, starting from rest. Af- 
ter a certain amount of acceleration, the vortex may be- 
come unstable, so that all the additional energy goes into 
producing spin waves. Alternatively, the vortex may re- 
main stable, and absorb all the energy. The energy goes 
into stretching the components apart until the vortex be- 
comes needle-shaped. The velocity is bounded in this 
case as well. E is proportional to the vortex-length D (x p 
and a linear dispersion implies a finite velocity. 

Spin-wave dissipation seems to be the fate of a vor- 
tex in a condensate where the atomic interactions in 
the Hamiltonian are rotationally symmetric (of form 
+ f /3(V'^SV')^). This conclusion is based mainly 
on numerical solutions of the Gross-Pitaevskii equations 
for steadily moving vortices with /3 = —.3a,q = .5/i. 
(see Fig. [3]). The computer did not find solutions past 
Vc = .65 y/^, which is close to the spin wave speed. 

Experimental conditions can be adjusted so that the 
vortex does not radiate sound. The symmetry of the 
condensate has to be broken further, by displacing the 
traps for the three Sz states into parallel planes. They 
still have to overlap some to allow for interspecies conver- 
sion. This trap set-up maintains the distribution of the 
atoms among the spin states more rigidly [lo]. The com- 
ponent vortices can now separate arbitrarily far without 
any instability. 

Displaced clouds also allow one to use the component- 
vortex picture more rigorously [l^. When \f3\ <C amin 
{oimin is the Smallest eigenvalue of ami 7712)1 the phase- 
only approximation j^/iml = const, applies. One can then 
derive the p{D) relation Eq. ((T^) and show that the speed 

of an infinitely stretched vortex is \J 

To summarize, each charge-vortex has only two spa- 
tial degrees of freedom because the lift force overcomes 
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FIG. 3; Computer solution for vortex cores moving at speeds 
w = 0, .45^/^, with rotationally symmetric interactions. The 
darker regions have the higher energy densities. The arrows 
indicate where ipi and ip-i vanish. Note that the energy 
density extends past these zeros, where the phase mismatch 
ends. This energy must come from varying magnetization, 
which eventually causes spin wave emission. 



the inertia. In contrast spin vortices in an unmagnetized 
condensate behave like classical particles because they 
are made up of oppositely rotating vortices whose total 
Magnus force cancels. The internal stretching between 
these components give rise to the mass. Spin vortices 
have four degrees of freedom: the center of mass coor- 
dinates and the momenta which are proportional to the 
distortion of the vortex core. The vortices behave like 
classical particles at low speeds, but betray their com- 
posite origin when accelerated sufficiently. 
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